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Abstract 

Among tlie firm predictions of softiy broken supersymmetry is the identity of gauge 
couplings and the corresponding Yukawa couphngs between gauginos, sfermions and 
fermions. In the event that a SUSY-hke spectrum of new particles is discovered at fu- 
ture colliders, a key follow-up will be to test these relations experimentally. In detailed 
studies it has been found that the SUSY- Yukawa couplings of the electroweak sector 
can be studied with great precision at the ILC, but a similar analysis for the Yukawa 
coupling of the SUSY-QCD sector is far more challenging. Here a first phenomenolog- 
ical study for determining this coupling is presented, using a method which combines 
information from LHC and ILC. 



1 Introduction 



With more than 1 fb~^ of luminosity on tape at the Tevatron experiments, and with the LHC 
and ILC on the horizon, TeV scale physics is fast entering the realm of experimental study. 
Among the most interesting and comprehensively studied possibilities for observable new 
physics is supersymmetry (SUSY). Interesting in its own right by virtue of being the largest 
possible space-time symmetry [1], it was chiefly with a number of additional theoretical 
successes in the early eighties that supersymmetry gained widespread acceptance, among 
these successful gauge coupling unification, radiative breaking of electroweak symmetry, and 
a natural candidate for Dark Matter. 

Stated briefly, supersymmetry promotes all the flelds of the Standard Model (SM) to su- 
perfields, all its multiplets to supermultiplets. Each such multiplet contains one boson and 
one fermion, which are related to each other by supersymmetry. Hence not only should all 
the SM particles (including at least one more Higgs doublet) have 'sparticle' superpartncrs, 
with spins differing by 1/2, but also the interactions of these superpartners should be fixed 
by supersymmetry. Even when considering broken supersymmetry (the only phenomeno- 
logically viable option), many of these relations persist, or are only modified to a tractable 
degree. Hence, part of the effort of testing the hypothesis, should suitable particles be dis- 
covered, will be to check whether their interactions do obey the supersymmetric relations or 
not. It is the feasibility of a particular such study we shall address here. 

Among the most intuitively clear relations is that the supersymmetric partner of a gauge 
boson, a gaugino, must couple with the same strength to gauge charge as its partner does. 
That is, the Yukawa coupling g between a gaugino interacting with a fermion and a sfermion 
must be identical to the corresponding SM gauge coupling g between a gauge boson and 
two (s)fermions. More specifically, denoting an SM fermion (gauge boson) by / {V), and its 
sfermion (gaugino) superpartner by / (V^), 

g = g{Vff) = g{Vff), 9 ^ g{Vff) , (1) 

with 

g-9 (2) 

required by supersymmetry. Indeed this relation is vital for the cancellation of quadratic 
divergencies in the radiative corrections to the Higgs mass. 

In what follows, we now constrain our attention to the Minimal Supersymmetric Exten- 
sion of the Standard Model (MSSM), with CP and i?-parity conserved. Hence the Lightest 
Supersymmetric Particle (LSP), usually the lightest neutralino, is stable, and supersymmet- 
ric particles can only be produced in pairs. 

If the supersymmetric particles are light enough to be produced at the next generation of 
colliders, it has been shown [2] that the supersymmetric Yukawa couplings in the electroweak 
sector can be precisely tested at a high-energy e+e^ collider. The t-channel contributions 
to neutralino/chargino [3] and slepton production [2,4,5] in e'^e~ and e~e~ collisions are 
directly dependent on the U(l) and SU(2) Yukawa couplings. Prom measurements of total 
and differential cross sections, these couplings can thus be extracted with a precision at the 
per-cent level [2,4,5]. 
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However, the analysis of the supersymmetric Yukawa couphng Qs in the SU(3) QCD 
sector is much more difficult. At the ILC this interaction can be studied in the process 
e+e~ (10*9-1 QQg [6]. While this approach is in principle straightforward, it suffers both 
from a very small cross section, at most about 1 fb for squark and gluino masses of a few 
hundred GcV [6], and also from the daunting challenge of selecting the signal from much 
larger ti and squark backgrounds. At the LHC on the other hand, squarks and gluinos 
with masses below 2-3 TeV are copiously produced, and their pair production cross sections 
depend directly on the supersymmetric Yukawa coupling Qs- However, measurements of total 
cross sections are exceedingly difficult in this environment, with typically only one or two 
specific decay channels of the squarks and gluinos experimentally accessible [7]. 

Here we consider a combination — the relevant branching ratios are to be determined at 
an e~^e~ collider and combined with exclusive cross section measurements in selected channels 
at the LHC. Together, this information determines the total squark or gluino production cross 
section, from which a value for the SUSY-QCD Yukawa coupling Qs can be extracted. 

In this publication, a ffist analysis is presented to explore whether such a coherent 
LHC/ILC analysis can provide a useful determination of the SU(3) Yukawa coupling. Since 
it is not clear a priori that the goal is achievable at all, we constrain our attention to a 
rather optimistic scenario in this study, with reasonably low masses, large cross sections, 
and a fairly large squark-gluino mass splitting. 

In the next section, the production of squarks at the LHC is studied and the basic 
phenomenology laid out. The relevant branching ratios of the squarks and how they could 
be determined from squark pair production at e+e" colliders is the topic of section 3. We 
round off by presenting the combined result (in the particular scenario studied here) and 
give a few concluding remarks. 

2 Squark production at the LHC 

2.1 Phenomenology and Strategy — LHC 

In pp collisions, squark and gluino production occurs in a variety of combinations, often with 
both s- and t-channel diagrams as well as several different partonic initial states contributing 
to the same final state. Fig. 1 contains a condensed summary of these, with explicit illus- 
tration of which vertices go with the QCD gauge (dots) and Yukawa (squares) couplings, 
respectively. Naively, there are several processes from which one might attempt to extract 
the Yukawa coupling. However, final states which receive contributions from both gauge and 
Yukawa couplings make the analysis more complicated and model-dependent. To isolate the 
Yukawa coupling, we note that the production of same-sign squarks (e.g. ulUl, uru^i, ...) 
proceeds only through the diagram shown in the lower left corner of Fig. 1, if one neglects 
the much smaller electroweak contributions. Hence, to excellent approximation, this pro- 
cess depends solely on the supersymmetric Yukawa coupling gs, and to the extent it can be 
isolated from the background, a clean determination of Qs should be possible. Note that, 
since the final state flavours are locked to the initial state ones, in pp collisions this process 
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Figure 1: Some examples for Feynman diagrams for partonic squark and gluino production 
in hadron collisions. Dots indicate the gauge coupling g^, while squares stand for the Yukawa 
coupling Qs. 



dominantly produces u and d squarks, with smaller admixtures of sea-flavoured squarks, in 
direct proportion to the quark content of the proton at the relevant x and values. 

Due to the flavour locking, only the first two generations of squarks are thus relevant, 
for which mixing effects arc small and we can take mass and current eigenstates to be 
identical to good approximation. That is, the heavier q mass eigenstate is pure gL (weak 
isospin doublet), and the lighter one pure gR (weak isospin singlet). Nominally, the lighter 
one would be the better target for a high-statistics study, simply due to phase space, but 
since it doesn't couple to weak interactions, it decays almost exclusively via the hypercharge 
coupling to a same-flavour quark and the LSP. Since charge tagging for light-flavour jets is 
exceedingly difficult, this decay mode effectively obscures the fact that we had same-flavour 
squarks to begin with. Moreover, since it only contains a jet and missing energy, the mode 
would be extremely challenging to separate from the background. The only feasible avenue 
thus appears to be to use flavour/charge tagging modes of the heavier mass eigenstates, the 

For ^L, the charge of the squark can be tagged through a chargino decay chain, 

ul^ dxt ^ dl'^ I'lXi: dj^^uxT ^uTPiXi: (3) 

dxT ^ dr PiXi, dl^uxt ^ul+uiXi, (4) 

and similarly for sl and cl. For a given squark flavor, the sign of the final-state lepton is 
related to the charge of the (anti-)squark. The production of same-sign squarks through 
the diagram in the lower left corner of Fig. 1 with this decay channel will therefore lead 
to same-sign leptons in the final state, while other direct squark production processes will 
tend to produce opposite-sign leptons in the final state. At this level, the signal is thus 
characterized by two same-sign leptons, two hard jets and missing transverse energy in the 
final state. 
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A very problematic background can come from gluino pair and mixed gluino-squark 
production if rrig > rriq^. In this case, gluinos can decay into quarks and squarks, g qqh, 
generating a component of two same-sign squarks plus additional jets which has a much 
more complicated (see Fig. 1) dependence on both the QCD gauge coupling and the Yukawa 
coupling we are interested in. This background becomes particularly challenging if the mass 
difference rrig — mq^ is small, since then the additional jets from the gluino decay will be soft. 
Due to the large overall mass scale in the process, the bulk of the cross section will contain 
several extraneous QCD jets, e.g. from initial state radiation, which can be significantly 
harder than this [8], and hence the separation of squark production into direct and gluino- 
induced samples will be extremely non-trivial. 

In the following we will only consider a scenario where the rrig — rrig^ mass difference is 
sufficiently large to allow a veto on additional jets for gluino background reduction. Although 
this will also reduce the signal somewhat, there are indications [8] that the particular signal 
process considered here, being dominated by valence quarks in the initial state, is associated 
with less extra radiation than other SUSY processes, which generally (in pp) have one or 
more gluons and/or sea quarks as initiators. 

The most important backgrounds from Standard Model sources are W^W^jj (same- 
sign VF's), where j is a light-flavour jet, and semi-leptonic tt, with the second lepton coming 
from the decay of a bottom quark. Due to the large total tt cross section, this can result 
in a sizable background. Below we include both these sources in our estimates, though the 
level of sophistication is obviously not anywhere near what would be required for a real 
experimental analysis. 



2.2 Numerical Results - LHC 

For the numerical analysis we use the scenario given in the appendix. This scenario is 
similar to the Snowmass point SPSla [9], with the exception that the gluino mass is raised 
to 700 GeV. Since tan P = 10, the lightest scalar tau eigenstate has a sizeable left-handed 
component, and since the other left-handed sleptons are too heavy, the chargino mainly 
decays into scalar taus, which subsequently decay into taus. To trace the charge explicitly, 
we here restrict our attention only to the leptonic tau branching fraction. The decay chain 
for the signal process is then 

65% J 100% , „ 35% 

ul ^ dxt >dT^Vrx\ >dt + :^, 

61% _ 100% „ 35% 

di, > ux-^ > UT frXi +-^, i = e, fi. (5) 

The numbers above the arrows indicate branching fractions and ^ stands for missing energy. 

Both signal and top and gluino backgrounds were simulated with Pythia 6.326 [10], 
i.e. with the hard 2 —>■ 2 scattering process calculated at leading order, dressed up with 
sequential resonance decays, parton showers (virtuality-ordered 'power' showers), underlying 
event ('Tunc A' [11]), string hadronisation, and hadron decays (the tt'' was set stable to 
optimize efficiency). The WWjj background was generated with MadEvent [12], i.e. at 
leading order with stable outgoing W bosons and no fragmentation. 
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The cross sections for squark and gluino production were normalized by the iC-factors 
obtained with Prospino 2.0 [13], which includes next-to- leading order QCD corrections. 
Top production was likewise normahzed to a total inclusive LHC cross section of 800 pb, 
while for the W^W^jj background only leading order results are available, translating to a 
comparatively larger theoretical uncertainty on this contribution. We further used Pythia 
to estimate the possible contamination from VV {ZZ, W^Z, and production, and 

found it to be small. 

Finally, a mock detector was set up, based on a toy calorimeter spanning the pseudora- 
pidity range \r)\ < 5 with a resolution of 0.1 x 0.1 in x space and a primitive UAl-like 
cone jet algorithm, with a cone size of AR = 0.4. Muons and electrons were reconstructed 
inside \ri\ < 2.5 and an isolation criterion was further imposed, requiring both less than 10 
GeV of additional energy deposited in a cone of size AR = 0.2 around the lepton and also 
no reconstructed jets with ptj > 25 GeV closer than AR — 0.4 around the lepton. These 
criteria duphcate the default settings of the ATLFAST simulation package [14] . 

Based on the general characteristics of the squark signal, we chose a set of preselection 
cuts, to broadly define the signal region: 

• at least 100 GeV of missing energy. 

• at least 2 jets with pxj > 100 GeV. 

• Exactly two isolated leptons i — e, /i with pt,^ > 7 GeV. 

As evident from Tab 1, after preselection most backgrounds are still much larger than 
the same-sign squark signal. The backgrounds can be further reduced by making use of the 
following characteristics. 

2.2.1 Cuts — Set A 

The signal is flavour-locked to the initial state and hence contains almost no heavy-flavour 
squarks. On the other hand, both tt and gg backgrounds contain large b fractions (it is a gen- 
eral feature of models with universal high-scale squark masses that gluinos decay dominantly 
into third-generation squarks). 

The number of b jets (inside \ri\ < 5 and with pt > 25 GeV) which would be identified 
with a perfect b tagger is plotted in Fig. 2a. In reality, b tagging is a trade-off between 
efficiency e and mistagging rate D. Following the ATLAS Physics TDR [14] we investigated 
the performance of three different b taggers; one. Fig. 2b, with a large efficiency e = 90% 
making it efficient against the background, but also high mistagging D = 25% rate, making 
it expensive on the signal, a second. Fig. 2c, with e = 80% and D — 10%, and finally a third 
with e = 60% and D = 2%. 

For our purposes, the most important goal is to obtain as pure a signal as possible, and 
hence we choose the Ob bin of the highest efficiency tagger. Fig. 2b, for this analysis. For 
completeness, note that restricting the b jet algorithm to the central detector region \r)\ < 2.5 
would not greatly affect our results; only a very small fraction of the true b jets lie outside 
this region. 
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Figure 2: Distribution of number ofb tags for the 4 different taggers discussed in tlie text, 
(a) being an idealized tagger and (b)-(d) ATLAS-inspired ones. 




Figure 3: (a) Distribution of the missing transverse energy after cuts on b-tagging. The 
bins below 100 GeV are cut out by the preselection, (b) Distribution of the transverse 
momentum of the third jet after the cut on ^. 
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Cross Sections Signal Backgrounds 



Y^a-udsc'^ (fb) 




Sum 


tt 








99 


'"P -i- 1 

iotal 


2100 




8x10 




7000 


1350 


3200 


Preselection 


49.2 


384.6 


177.7 




136.4 


23.2 


47.3 


b-veto 


17.1 


31.4 


13.0 




10.3 


7.1 


1.0 


^ > 150 GeV 


15.1 


22.2 


6.1 




9.0 


6.2 


0.9 


A) pTj3 < 75 GeV 


10.1 


9.5 


3.6 


N/A 


2.3 


3.5 


0.1 


Ptji > 200 GeV 


8.9 


<7.6 


1.7 


<0.7 


2.0 


3.2 


0.1 


B) ptj3 < 50 GcV 


7.8 


5.9 


2.4 


N/A 


1.0 


2.5 


0.03 


PT,n > 200 GeV 


7.0 


<4.9 


1.0 


<0.7 


0.8 


2.3 


0.03 


iQe = +1) 


4.8 


<3.0 


0.6 


<0.5 


0.5 


1.4 


0.01 



Table 1: Total and progressively reduced cross sections estimated from the cuts deGned 
in the text, for the signal, summed background, and individual backgrounds, respectively. 
Two different sets of cuts have been considered, A and B, the latter incorporating a more 
severe 3rd jet veto. For illustration, the effect of including a further cut on the lepton sign 
is included for set B, though note that this last cut is not actually made in the analysis. 

Including now the b jet veto, Fig. 3 shows the distribution of missing energy for the 
remaining events. While the contribution from top still dominates the background it can be 
effectively suppressed by increasing the ^ cut. Here, we choose ^ > 150 GeV. 

At this point, the gluino-related backgrounds dominate. A veto on hard additional jets 
is essentially the only way of obtaining significant further rejection power against these. 
Fig. 3 (b) shows the distribution of the transverse momentum of the third-hardest jet for 
the events remaining in the analysis. By rejecting all events with pxjs > 75 GeV, the ratio 
of the signal to gluino background is markedly improved. 

The largest remaining background is now top production again (closely followed by qhQl), 
whose two hard jets tend to be softer than those of the signal, due to the large squark- 
chargino mass difference. This is illustrated in Fig. 4 (a), where the transverse momentum 
of the hardest jet in the remaining events is plotted. 

Increasing the transverse momentum cut on the first jet to Pt,\i > 200 GeV wc obtain a 
final signal cross section of 8.9 fb, over a background of less than 7.6 fb, with the effect of 
progressive cuts on signal and background normalisation and composition given in Tab. 1. 
The signal-to-background ratio is 1.2, sufficient to allow a meaningful measurement. With an 
integrated luminosity of 100 fb~^, the statistical error on the same-sign squark cross section 
is 4.5%. 

The reason we only quote an upper estimate of the background in Tab. 1 is that the 
WWjj sample was not improved by parton showering, hence we have no way of estimating 
the effects of the 3rd jet veto on this component. 
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Figure 4: Distribution of the transverse momentum of the hardest jet after the third jet 
veto for (a) set A, and (b) set B. The bins below 100 GeV are cut out by the preselection. 



2.2.2 Harder Cuts — Set B 

It is interesting to check if the situation can be improved by choosing stricter selection cuts. 
The gluino background can be further reduced by tightening the cuts on the third jet to 
Ptjs < 50 GeV. In addition, the signal is dominated by up-squarks, since two out of three 
valence quarks of the proton are up-quarks. The qq* background on the other hand needs 
at least one sea-quark in the initial state and therefore results in up- and down-squarks 
in more equal proportion. Thus the requirement of positively charged leptons in the final 
state will increase the ratio of same-sign squarks over opposite-sign squarks. With these 
stronger cuts the cross section values on the last line of Tab. 1 are obtained. Indeed the 
signal-to-background ratio is improved to 1.6, but at the cost of reduced overall statistics. 
We choose a compromise; omitting the final cut on lepton sign results m S/B = 1.45, with 
acceptable statistics, yielding a signal cross section of 7 fb, over a total background of less 
than 4.9 fb. With a 100 fb~^, the statistical error on the signal is 4.9%. This error is slighty 
larger than the result of the cuts Set A, but this is offset by lower background systematics 
(see section 4). 

Noting that the two analyses are thus identical apart from the 3rd jet veto, one might 
ask why it is important to show the results of both. The reason, slightly touched on above 
and discussed in more detail in [8], is that 50 GeV jets are quite soft objects at the LHC, 
especially as compared to the signal scales we consider here, of order 600-700 GeV. Since 
we cannot be entirely sure that the parton shower resummation employed here is sufficiently 
accurate in this region, an explicit variation by a 25 GeV shift in the 3rd jet veto represents 
one way of verifying whether the analysis and thereby our conclusions are stable against 
variations of this cut. 



8 



3 Squark decays at an e+e collider 
3.1 Phenomenology and Strategy — ILC 

In order to interpret the measured rates at the LHC in terms of the total squark production 
cross section, the individual branching ratios in the decay chain eq. (5) need to be established. 
We now turn to a discussion of how each of these could be extracted from measurements at 
a high-energy e^e~ linear collider. 

The chargino properties are the most straightforward to determine, being easily accessible 
via the chargino pair production process, e+e~ XiXi- Due to the large rate for that 
process, the expected error on cross section measurements is about 1% at a 500 GeV collider 
[7]^. Since the charginos xf are among the hghtest supersymmetric particles, the only open 
decay channels are x^ I'^^iXi^ ^ — ^tI^^^^ which can all be easily separated from the 
background [7]. We therefore assume that the chargino branching ratios are known with 1% 
error: BR{xt ^ r+UrXi) = (100 ± 1)%. 

The squarks offer substantially more resistance. In the given scenario the left-squarks 
are in fact shghtly too heavy to be accesible at a linear coUider with 1 TeV center-of-mass 
energy. Rather than attempting to find a more amenable parameter point, we here consider 
the hypothetical case of an e~^e~ collider with a center-of-mass energy of about 1.5 TeV. 
Obviously, we do not propose to address the technical challenges of such a machine, but 
justify our choice simply by the proof-of-concept nature of our study. We also assume that 
both beams can be polarized, with polarization degrees of 80% for electrons, and 50% for 
positrons. 

Having kept the squarks heavy, we must now deal with the fact that they can decay into 
the whole spectrum of charginos and neutralinos. The heavier charginos and neutralinos have 
large branching ratios into gauge bosons, which distinguishes them from the lighter states. 
The following characteristic decay modes are chosen to distinguish the various charginos and 
neutralinos: 



xt ^^^rXl, 



Xi 



xl^TTxl 
xlA^W^xt^W^r^^rXl 



BR=100% 
BR=24% 
BR=100% 
BR=59%,52% 



(6) 
(7) 
(8) 
(9) 



A distinction of the and X4 states is not necessary since for the purpose of this study only 
the branching fraction into the lightest chargino needs to be known as an absolute number. 
The tau leptons in the final state can be identified in their hadronic decay mode with roughly 
80% tagging efficiency. Since about 65% of tau decays are hadronic, this amounts to a total 
tau tagging efficiency of about 50%, consistent with the findings e.g. of LEP2 studies (see 
for instance Ref. [15]). 

^Notc that the SPSla scenario studied in Ref. [7] is identical to the scenario in the appendix in the 

electroweak sector. 
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3.2 Numerical Results - ILC 



For this work, Monte-Carlo samples for squark pair production in the different squark decay 
channels have been generated at the parton level with the tools of Ref. [5, 16]. Also the most 
relevant backgrounds have been simulated, stemming from double and triple gauge boson 
production as well as tt production. It is assumed that an integrated luminosity of 500 fb~^ 
is spent for a polarization combination P{e^)/P{e') — +50%/— 80%, which enhances the 
production cross section both for Ml and c^l production. Here =F indicates left /right-handed 
polarization. The branching fractions are obtained from measuring the cross sections of all 
accessible decay modes of the squarks and identifying the fraction of decays into one specific 
decay mode out of these. 

Since the squarks are produced in charge-conjugated pairs, it is a priori difficult to dis- 
tinguish up- and down-squarks in the final state. However, assuming universality between 
the first two generations, a separation between up- and down-type squarks can be obtained 
through charm tagging. According to Ref. [17], a c-tagging efficiency of 40% is achievable 
for a purity of 90%. By combining the different decay channels in eqs. (6-9), the following 
final state signatures are identified as interesting: jj{nT)^ with n e {1,2,3,4}, cc{nT)^ 
with n G {2, 4}, jjrr^Z/W)^, ccttZ^, where j indicates an untagged jet, c a tagged charm 
jet, and Z/W a hadronically decaying gauge boson where the invariant mass of the two jets 
combines to the given gauge boson mass. For charged current squark decays, charm tagging 
does not provide any additional information, since a cs-|-leptons final state can arise both 
from up- and down-type squarks. Therefore, only signatures with two tagged charm jets are 
included in the list above. Finally, since several squark decay channels can contribute to 
most of the final states above, one has to solve a linear equation system in order to derive 
the individual contributions. 

Before presenting results for this procedure, there is an additional complication due to the 
subsequent decays of the neutralinos and charginos after the squark decays. The branching 
fractions for the light chargino xf neutralino X2 can be deduced from cross section 
measurements at a low-energy run with ^/s = 500 GeV, as mentioned above. Only a few 
decay channels are kinematically open, and the environment is sufficiently clean so that all 
decay channels are easily visible with a small error. 

The situation is more tedious for the heavier charginos and neutralino, which have many 
different decay channels open. Rather than attempting to measure individual cross sections 
for all these decay modes, here instead a different method is pursued. Near production 
threshold of two supersymmetric particles XY, the energy of a particle from a two-body 
decay mode of e.g. X is sharply defined. The signal of such a monoenergetic particle would 
clearly stand out over the background, which is smooth in the particle energy, and can thus 
be used to identify a particular production cross section irrespective of the decay mode of 
the other sparticle Y. 

The branching ratios of the heavy neutralinos X3 4 can be studied in associated production 
with X2: e+e" — > X2X3,4- When approaching the production threshold, the decay of the X2 
is characterized by a monoenergetic r lepton over a smooth background. The production 
threshold for X2X3 ^i^s at ^/s — m^o + m^o, 535 GeV in the present case. A few GeV 
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uxl 
dxt 
dxi 



0.9 ± 0.5 % c^L 
29.0 ± 3.0 % 
<1% 

< 1 % 




1.9 ±0.8 % 
28.3 ± 4.8 % 



1.9 ±0.8 % 
63.9 ± 5.2 % 
4.0 ± 1.4 % 



< 0.2 % 



67.7 ±3.2 % 
1.4 ±0.7 % 



Table 2: Estimated error for squark branching determination from measurements at an 
e+e~ collider with ^/s —1.5 TeV. 

above this threshold, at ^/s = 540 GeV, the production cross section is already sizeable, 
a23{y/s = 540 GeV) = 15.8 fb with polarized beams P(e+)/P(e-) = ±50%/-80%. The 
cross section for X2X4 is unfortunately too small to allow a meaningful measurement. An 
alternative is to obtain the branching fractions of Xa the XsXa threshold, by tagging 

the X3 decays into a monoenergetic 1^-boson. Of course for this purpose only hadronic W 
modes can be used. A few GeV above the kinematical threshold m^o + m^o = 737 GeV, 



the production cross section is relatively large, a34{^/s = 745 GeV) = 29.3 fb, assuming 
polarized beams as above. 

For a full feasibility demonstration of this technique, a realistic detector simulation would 
be necessary. Here the achievable precision for the neutralino and chargino branching ratios 
is estimated by simulating the parton-level production near threshold. Assuming 50 fb~^ 
luminosity for each of the threshold measurements, the following statistical errors are found: 



Together with this information, the squark cross section measurements at y/s — 1.5 TeV 
can be interpreted in terms of squark branching ratios. Solving the linear equations system 
that connects final state signatures with production processes with a x^ scan, the estimated 
errors for the squark branching ratios are listed in Tab. 2. For the purpose of this work, the 
interesting information are the branching ratios of and c^l into the light chargino xf- 



Based on the results from the simulations for squark production at the LHC and the ILC 

presented in the previous sections, one can now derive an estimate for the precision for the 
determination of the supersymmetric strong Yukawa coupling Qs. Since the production of 
same-sign squarks at the LHC with gluino t-channcl exchange is proportional to cfgLgL] oc gf, 
the error for the coupling constant Qs is roughly one quarter of the error of the total same- 
sign squark production cross section at the LHC. The statistical uncertainty is combined 
with the most important systematic error sources in Tab. 3. The numbers are based on the 
choice (B) for the LHC signal selection in Tab. 1. We now turn to a discussion of the various 
unvertainties and caveats associated with these numbers. 



BR(x^ ^ W^xf) = (59 ± 6.5)%, BR(x^ ^ W±xf ) = (52 ± 2.5)%, 
BR(x^^^xf) = (24±1.3)%. 



(10) 



4 Combination - LHC and ILC 
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9s/ gs 


LHC signal statistics 


4.9% 


1.3% 


SUSY-QCD Yukawa coupling 






m qi^g background 


2.4% 


0.6% 


PDF uncertainty 


10% 


2.4% 


NNLO corrections 


8% 


2.0% 


Squark mass Am^^ = 9 GeV 


6% 


1.5% 




8.2% 


2.0% 




17.3% 


4.1% 



Table 3: Combination of statistical and systematic errors for the determination of the 
absolute same-sign squark production cross section at the LHC and the derivation of the 
strong SUSY-Yukawa coupling. The cuts set B were used for the LHC 



The remaining background after selection cuts at the LHC introduces a systematic un- 
certainty since the supersymmetric background from gluino production also depends on Qg. 
Including the variation of gs in the dominant squark-gluino background introduces an un- 
certainty of 5% on the signal. The background from opposite-sign squark production is also 
sizeable, but depends in the same way on as the signal, a (x and thus does not introduce 
an additional systematic effect. 

In order to determine absolute cross sections, the proton parton distribution functions 
(PDFs) in the relevant x and ranges must be known accurately. The production of same- 
sign squarks proceeds from same-sign quarks in the initial state and is therefore dominated by 
the contributions from valence quarks. The value of valence quark PDFs at high scattering 
energies can be computed reliably from perturbative parton evolution and can be checked 
against measurements of vector boson production. We tested this assumption by comparing 
results for different CTEQ PDFs releases (CTEQ6M, CTEQ6D, CTEQ5M1, CTEQ5L and 
CTEQ6L1) [18] and find a maximal variation of 7% in the signal and 4% in the background. 
Thus we assign a total error of 10% for the parton distributions. 

Squark production receives large higher order QCD corrections. The next-to-leading 
order corrections are known [13] and can be included in the analysis. The uncertainty of 
the missing 0(0"^) contributions arc estimated by varying the renormalization scale of the 
0{as) corrected cross section between mq^^/2 < Q < 2mg^^, leading to an error of 8%. 
Furthermore the cross section depends on the values of the squark masses. According to 
Ref . [7] , the left-chiral first generation squark masses can be determined with an error better 
than Am^L = 9 GeV, leading to an uncertainty of about 6% for the production cross section. 
Finally, the expected error for the determination of the squark branching ratios at the linear 
collider must be included. 

Combining all error sources in quadrature, it is found that the same-sign squark produc- 
tion cross section can be determined with an error of 17.3%, translating to an error of 4.1% 
for the supersymmetric QCD Yukawa coupling gg. That is, assuming the steps outlined here 
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can be carried out in more or less the fashion described, it should be possible to test the 
super symmetry identity between gauge and Yukawa couplings to the level of a few percent. 

It should be noted though that this result is highly scenario dependent. For different su- 
persymmetric scenarios, with larger squark masses, smaller gluino masses or different squark 
decay modes, the analysis could be much more difficult or even impossible. It is also interest- 
ing to study gluino production not only as a background, but as a signal process for studying 
the SUSY-QCD Yukawa coupling, although this introduces additional complications since 
gluino production always depends both the the gauge and Yukawa couplings. Nevertheless 
we hope this study shows that the determination of the supersymmetric QCD Yukawa cou- 
pling with competetive precision is not out of reach, and that it is a tantalizing target for a 
complex all-embracing combination of LHC and lepton collider data. 

5 Conclusion 

We have described a phenomenological method for testing the identity between the QCD 
gauge coupling and the corresponding squark-gluino-quark Yukawa coupling that arises in 
supersymmetric theories. Noting that same-sign squark production at a hadron collider is 
both a very clean signal and also directly proportional to the fourth power of the QCD 
Yukawa coupling, we propose to measure the exclusive cross section for same-sign squark 
production at the LHC in a particular decay channel. To convert this to a measurement 
of the Yukawa coupling, branching fraction determinations at a future lepton collider are 
(de-) convoluted with the LHC measurement to yield a combined determination of the total 
inclusive cross section for same-sign squark production. We find that this cross section can 
be determined to roughly 20% precision, in the scenario studied here, translating to roughly 
5% precision on the Yukawa coupling. This is certainly encouraging and, we hope, sufficient 
in itself to spark off more in-depth studies. Obviously, due to the complexity and breadth 
of the proposed analysis, we were forced to employ some shortcuts here, as described in 
detail in the text above. It would be highly interesting to establish how these estimates 
would compare to both a more detailed (in the sense of detector capabilities and systematic 
uncertainties) and comprehensive (in the sense of exploring a wider range of parameter space) 
analysis. 

Appendix 

Here the reference scenario is listed that is used in this analysis. It is coincides with the 
Snowmass point SPSla [9], but has a larger gluino mass of 700 GeV. The scenario is defined 
at the weak scale through the supersymmetry breaking parameters. The parameters relevant 
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for this work are 

Ml = 99 GeV mL3 = 197 GeV mqi = 540 GeV 

M2 = 193 GeV mR3 = 136 GeV mui = 522 GeV 

M3 = 700 GeV Ar = -254 GeV mm = 520 GeV 

/X = 352 GeV tan p = 10. 

The tree-level masses are computed from this to 

= 537 GeV m^o = 96 GeV 

mj^ = 543 GeV m^o = 177 GeV 

= 133 GeV m,-,± = 176 GeV 



Xi 



= 700 GeV m^o ~ 360 GeV. 



;ii) 
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